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We calculate the finite temperature compressibility for two-dimensional semiconductor systems, 
monolayer graphene, and bilayer graphene within the Hartree-Fock approximation. We find that 
the calculated temperature dependent compressibility including exchange energy is non-monotonic. 
In 2D systems at low temperatures the inverse compressibility decreases first with increasing tem- 
perature, but after reaching a minimum it increases as temperature is raised further. At high 
enough temperatures the negative compressibility of low density systems induced by the exchange 
energy becomes positive due to the dominance of the finite temperature kinetic energy. The inverse 
compressibility in monolayer graphene is always positive and its temperature dependence appears 
to be reverse of the 2D semiconductor systems, i.e., it increases first with temperature and then 
decreases at high temperatures. The inverse compressibility of bilayer graphene shows the same 
non-monotonic behavior as ordinary 2D systems, but at high temperatures it approaches a constant 
which is smaller than the value of the non-interacting bilayer graphene. We find the leading order 
temperature correction to the compressibility within Hartree-Fock approximation to be In T at 
low temperatures for all three systems. 



I. INTRODUCTION 

In this article, we provide a detailed theory for the 
temperature-dependent electronic compressibility {K) of 
an interacting two dimensional (2D) electron (or hole) 
system within the Hartree-Fock theory. We consider 
three distinct 2D systems of active current interest in con- 
densed matter physics: 2D semiconductor systems (e.g., 
quantum wells, heterostructures, inversion layers, etc.); 
monolayer graphene; bilayer graphene. For the purpose 
of comparison (and also for the sake of completeness), 
we also provide as an appendix the corresponding finite 
temperature Hartree-Fock compressibility for a standard 
three dimensional electron gas (3DEG) since this result 
does not appear to be available in the theoretical liter- 
ature in spite of the sixty-year long history of studying 
many-body interaction effects in 3DEG. 

Ever since the pioneering measurement of the 2D com- 
pressibility by Eisenstein et al-r'^ it has been extensively 
studied in 2D systems because the compressibility pro- 
vides fundamental insight into quantum ground state 
properties which are not readily obtained from trans- 
port measurements. Some experimental studies"^"^ of 
the compressibility have shed light on understanding the 
2D metal-insulator transition (MIT)^. The inverse com- 
pressibility {K~^) is positive when the kinetic energy 
dominates over interactions at high densities. As the 
2D density is reduced changes sign and becomes 

negative due to the increase of exchange energy associ- 
ated with electron-electron interactionii^. The negative 

reaches a minimum value at a certain low density 
n, and then increases dramatically with further decreas- 
ing n as the 2D MIT sets in. The minimum point in 

has occasionally been loosely identified as the criti- 
cal density for the 2D MIT and is closely related to the 
transition from a homogeneous system to an inhomoge- 
neous nonlinear screening regimej^i^i where the system is 



dominated by the formation of electron (or hole) pud- 
dles. Whether this observed low-density behavior of K 
is a cause for or an effect of 2D MIT is unclear. The 
low-density compressibility of 2D electronic (in this pa- 
per we use the term "electron" to imply either electron or 
hole as should be obvious from the context) systems has 
been of interest for almost 20 years now with the early 
experimentsi"— studying semiconductor-based 2D sys- 
tems extensively and very recent experiments studying 
monolayer and bilayer graphene. The reason for the fo- 
cus on the density dependence of compressibility is that, 
in general, quantum interaction effects increase monoton- 
ically with decreasing density in a Coulomb system (with 
monolayer graphene being an odd exception where the in- 
teraction parameter, the so-called graphene fine structure 
constant, is independent of carrier density^), as discussed 
above, leading eventually to the electronic compressibil- 
ity becoming negative at low, but easily accessible, den- 
sities due to exchange effects. 

It is curious, however, that in spite of this intense inter- 
est in the strongly interacting low-density 2D compress- 
ibility, there has been little research on the theoretical 
functional dependence of 2D compressibility on temper- 
ature. This is strange because the scale for the tem- 
perature dependence of any electronic property is the 
Fermi temperature Tp = Ep/ks, which invariably de- 
creases with decreasing density. Thus, temperature ef- 
fects on 2D compressibility become progressively more 
important even at a fixed temperature as the electron 
density goes down, since the important dimensionless 
temperature t = T/Tp increases with decreasing den- 
sity even if T is kept fixed. For example, 2D GaAs holes 
have Ti? < 1 K for a hole density of 10^° cm^^, which 
means that T = 1 K is effectively a high-temperature 
regime for 2D hole densities < 10^" cm^'^! Our work in 
this paper takes a first step in correcting this omission 
in the literature. Our results are important for interpret- 
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ing low-density 2D compressibility (even at relatively low 
temperatures) existing in the literature. 

The thermodynamic isothermal compressibility (or 
simply compressibility) of a system is defined as the 
change of pressure with volume, K = —V~^{dV/dP)N,T: 
where N is the particle number, T is the temperature, 
V is the system volume, P = —{dE/dV)N,T is the pres- 
sure, and E is the total energy of the system. Thus, 
K^^ = V [d"^ E / dV"^) . In the non-degenerate or classi- 
cal systems the velocity of thermodynamic sound, wq, of- 
ten provides a convenient experimental measure of com- 
pressibility, Vq = where m and n — N/V are the 
mass and the density of the system, respectively. In the 
quantum limit the compressibility can be obtained the- 
oretically by using the theorem of Seitz-, (i.e., by con- 
necting K with the change in chemical potential with 
density, d^/dn, where /i is the chemical potential and n 
is the carrier density of the system). Another method 
of evaluating K is through the compressibility sum rule. 
This is an exact relationship between the compressibil- 
ity and the long wavelength limit of the static dielectric 
function^*^. In 2D quantum systems the compressibility is 
often measured from the quantum capacitance^*^ which 
is proportional to dn/dfi oc K . A scanning single elec- 
tron transistor has also been used to directly measure the 
density dependent chemical potential^, and hence d/i / dn 
or dn/d^. We will mostly discuss in this work the be- 
havior of Kq/K, where ifg is tbe non-interacting com- 
pressibility, and Kq/K cx dfi/dn. We note that some- 
times dn/dfj, is referred to as the thermodynamic density 
of states since at T = and for noninteracting system 
dn/diJ, = Do, the noninteracting density of states. 

In this paper, we refer to K and as compressibil- 
ity and incompressibility (i.e. inverse compressibility) 
respectively, and also discuss results for dfi/dn, which 
is sometimes directly measured experimentally, remem- 
bering that this quantity is directly proportional to the 
incompressibility. A part of our theoretical motivation 
for exploring the temperature dependence of compress- 
ibility in 2D semiconductor systems arises from recent 
experiments conducted on GaAs-based 2D systemsii. 
In addition, there has been substantial experimental 
interest in the temperature dependence of graphene 
compressibility^^ We find that the calculated temper- 
ature dependent compressibility including exchange en- 
ergy is non-monotonic in temperature in both graphene 
and 2D semiconductor systems. However, their tem- 
perature behaviors appear to be reversed, i.e., the in- 
verse compressibility in graphene cx d^/dn increases 
first with increasing temperature and then decreases at 
high temperatures while the inverse compressibility in 
2D semiconductor systems decreases with temperature 
first and then increases at higher temperatures. We find 
that the leading order temperature correction to 2D 
within Hartree-Fock approximation (HFA) is In t at low 
temperatures for both systems {t = T/Tp). We also find 
that the BLG compressibility has very weak temperature 
dependence because the non-interacting kinetic energy is 



independent of temperature. 

To calculate the compressibility we use Seitz 's theorem 
which is given hy^ 

1 = n^d^iiT)/dn (1) 

where /i(T) is the finite temperature chemical potential 
and n the free carrier density. In general, the direct 
measurement of d^/dn provides information on the ther- 
modynamic many body renormalization (for example, 
Fermi velocity) arising from electron-electron interaction 
effectsi^ii^. Our goal here is to theoretically calculate 
the renormalized dfi/dn in graphene and 2DEG including 
exchange interaction effects, or equivalently in the HFA, 
which should be an excellent quantitative approximation 
for compressibility in 2D systems. Our calculated carrier 
density and temperature dependence of compressibility 
(or dfi/dn) can be directly compared to experimental 
measurements in 2D systems including graphene. 

The rest of the paper is organized as follows. In Sec. 
im we introduce the formalism that will be used in cal- 
culating compressibility in 2D semiconductor systems, 
which includes the model describing the Coulomb sys- 
tem, the procedure of calculation, and the numerical and 
analytical results of compressibility in 2D semiconductor 
systems. In Sees. [HT] and |IV] we present the theoreti- 
cal formalism and our analytical and numerical results of 
compressibility in MLG and BLG, respectively. Sec. IVl 
contains discussions and conclusions. 



II. COMPRESSIBILITY IN 2D 
SEMICONDUCTOR SYSTEMS 

In this section, we present the compressibility of 2D 
semiconductor systems (2DS) such as Si inversion lay- 
ers in MOSFETs, 2D GaAs heterostructures or quantum 
wells. In some of our numerical calculations we use the 
parameters corresponding to the n-GaAs or p-GaAs sys- 
tem, since these are the most-studied 2D system in the 
literature. When ever possible, we provide our results in 
dimensionless units for universal applicability. 

In the absence of interaction, the non-interacting fi- 
nite temperature chemical potential, iiq{T), is calculated 
through the conservation of the total electron density, 
i.e., 

n{T) = no = J deD{e)nFie), (2) 

where hq — n{T — 0), D{e) — gm/2'K is the den- 
sity of states with g being a degeneracy factor in- 
cluding both spin and possible valley degeneracy, and 
nF{f) = [e^'^■^~^^°> -\- \]~^ is the Fermi distribution func- 
tion. We take h = 1 throughout. By solving Eq. ([2]) 
self-consistently we have the 2D chemical potential at fi- 
nite temperatures 

^lo{T)/EF=tHe^" -I) (3) 
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where Ep = kp/2m = 2imQ/gm is the Fermi energy at 
T = and i = T/Tp with Tp '= Ep/kB- 

To include interaction effects in the chemical potential, 
we calculate the exchange self-energy contribution within 
HFA. We then have the interacting chemical potential as, 



= /iQ + T.{kp,Ep) 



(4) 



where Y,{kp,Ep) is the HF exchange self-energy calcu- 
lated at Fermi momentum kp. The HF exchange energy 
with the electron-electron Coulomb interaction is given 
byiS 



^cx(^) 



d^k' 



V{\k-k'\)np{ek' 



(5) 



with V{\k-k'\) = 2T:e'^ /K\k-k'\ is the 2D bare Coulomb 
interaction in the momentum space (k is the background 
lattice dielectric constant), ek = k^/2m is the single par- 
ticle energy (m is the effective mass of the particle). In 
the following, we calculate the compressibility of 2D semi- 
conductor systems for both zero temperature and finite 
temperatures by finding the total chemical potential in- 
cluding exchange effects as a function of density and tem- 
perature. 



A. Zero temperature compressibility 

At zero temperature {T = 0), the non- interacting 
part of the chemical potential is just the Fermi energy 
Ep. Then d^/dn becomes the inverse of the nonin- 
teracting single-particle density of states at the Fermi 
level: {dn/dfi) = Dq{Ep) — gm/2TT^ which is a density- 
independent constant for 2D. From Eq. ([T|) we have the 
zero temperature compressibility of non-interacting 2D 
systems, Kq, 



2 2 IP, 



1 



dn 



(6) 



The noninteracting compressibility, K^, is a positive 
quantity and is inversely proportional to the square of 
particle density, Kq cx n~^. 

For interacting 2D systems the exchange self-energy 
can be calculated at T = from Eq. ([Sj, 



Sox(fc) 



2e^ , f k \ 



(7) 



where E{x) is the the complete elliptic integral of the 
second kind. The total chemical potential within HFA 
can be calculated by putting k = kp, i.e.. 



fj. = fio + T,{kp) = Ep{l 



2V5 



(8) 



where g = g^gs {gs, gv are the spin and valley degen- 
eracy, respectively) and the dimensionless coupling con- 
stant Vs = 2e'^'m/ (y/gnkp) is the ratio of the average 



Coulomb potential energy U to the average T = ki- 
netic energy. The parameter as defined above is the 
dimensionless interaction parameter for a 2D electron liq- 
uid. We show most (except for Sec. HID) of our results in 
the dimensionless units of t = T/Tp and rs, which would 
apply to any 2D semiconductor system. Differentiating 
Eq. ([5]) with respect to n and making use of the rela- 
tion kp = Airn/g, we get the zero temperature inverse 
compressibility within HFA 



K 



Vgrs 



(9) 



The interacting inverse compressibility Kq / K monotoni- 
cally decreases as Tg increases (or density decreases) and 
changes its sign from positive to negative at a coupling 
strength rs = t^I y/g- This behavior has been observed 
in experimentsii^ and much discussed in the literature. 
Generally, negative compressibility leads to a thermody- 
namic instability of a system. However, the compress- 
ibility we have discussed in this paper applies only for 
the electronic part, i.e., the quantity is not actually the 
compressibility of the whole system including the positive 
background charge which is necessary for neutrality. The 
interaction between the electrons and its associated pos- 
itive neutralizing background has been ignored. In fact, 
a system with the negative electronic compressibility can 
be stabilized by the positive background, which gives rise 
to a positive compressibility for the whole system. The 
negative electronic compressibility of the 2DEG has been 
directly measured in experiments^"— which probe only 
the electronic part of the compressibility, obviously the 
total compressibility of such a system would still be pos- 
itive in order to maintain thermodynamic stability of the 
whole system. Results given in Eqs. dU-dl]) are of course 
well-knowi*^, and are given here only for the sake of com- 
pleteness and for the sake of comparison with our finite-T 
results discussed below. 



B. Finite temperature compressibility 

At finite temperatures the chemical potential with ex- 
change self-energy can be calculated from Eq. @ and the 
normalized chemical potential, jl — ^i/Ep is expressed as 
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dxKi [x] 



1/^2 



(10) 



where /io(T) — [ii^jEp is the non-interacting chemical 
potential given in Eq. t — T/Tp, and Ki{x) is the 
complete elliptic integral of the first kind. 

In the absence of exchange interaction, we obtain the 
non-interacting inverse compressibility at finite temper- 
atures 



1 



Km{T) 



dn Ko ' ^ ^ 
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where h{t) = ■^_^-i/t and K^i is the finite temperature 
non- interacting compressibihty. For T <C Tp, h{t) ^ 
1 + exp(-l/i). Thus we have K„,iT) - Ko{l - e^i/*). 
The leading order correction to Kni{T) at low tem- 
peratures is exponentially suppressed. For T ^ Tp, 
h{t) ^ t. Thus the noninteracting inverse compress- 



ibility increases linearly in the high temperature limit, 
K-l{T:^TF)^t/K„. 

Including exchange effects, the finite temperature in- 
verse compressibility Kq/K{T) can be obtained from 
Eq. (Uni) 



K{T) 



h{t) 



TT 
2TTt 



(lxKi(x) 



'lxKi(x) 



g(a;2-po)/t _|_ Y 

xix^ - h{t)) 



cosh 



-t- 



cosh' 



2 jl/x^-jio) 



(12) 



The full numerical results of the compressibility are given 
in Sec. Ill Cl We first consider the asymptotic behavior of 
the temperature dependence both at high temperatures 
(T > Tf) and at low temperatures (T < Tp)- 

In the low temperature limit {t <^ 1), the integrand in 
the temperature-dependent term is only appreciable near 
the divergent regime of the first elliptic function Ki{x) 
{x « 1). We could get the asymptotic form by expanding 
the elliptic integrals about that point, i.e.. 



Ki{x^\)^\n 



(1 



In 



-1 



(13) 

Replacing Ki{x) with Eq. ([T3|) in the integrand of 
Eq. PU)). the chemical potential in the low temperature 
limit T /Tp ^ 1 is given by 



2^5 



16. 



r,[l + -(2A + B-Aln— )] (14) 

TT O t 



where A = 7rVl2, B = f2[ln(47r) + 12C'(-1)] « 0.449 
and C'(~l) is the derivative of zeta function. Differenti- 
ating the asymptotic formula of the many-body chemical 
potential, we obtain the inverse compressibility for i <C 1, 



K{T) 



T^Tp 



hit) 



V~9. 



1 



0.136t2- — i2 



(15) 



Since ft,(i) ~ 1 for i ^ 1, the leading order correction 
to the inverse compressibility K~^{^) comes from the 
exchange energy term, i.e., t^lnf, which gives rise to 
the decrease of if o at low T with increasing temper- 
ature as shown in Fig. [Tfa). The asymptotic formula 
given in Eq. (jlSp agrees well with our numerical calcu- 
lation at low temperatures as shown in Fig. [lla). The 
low-temperature behavior of K^^{T) is dominated by 
the In t term which produces the shallow minimum in 
Kq/K as a function of i = T /Tp with the size of the mini- 
mum increasing with increasing (decreasing) (density) . 
This interesting low-temperature non-monotonicity in 
the temperature-dependent 2D compressibility is entirely 



an exchange effect. We emphasize, however, that this 
exchange- induced minimum is very shallow in Kq/ K{T), 
and never exceeds 10% of its T = value (often it is 
much less). 

The asymptotic behavior of high temperature com- 
pressibility is obtained by approximating the Fermi-Dirac 
distribution by the classical Boltzmann distribution. For 
t — > oo, the normalized chemical potential with exchange 
energy is given by 



il{T) = ilo - 



2Vt 



(16) 



The corresponding asymptotic formula for the high tem- 
perature inverse compressibility is given by 



K{T) 



t:s>Tf 



(17) 



For i ^ 1, the exchange energy contribution to the 
inverse compressibility decreases as i^^/^ while the ki- 
netic energy contribution increases lineally since h{t) ~ t. 
Thus, the kinetic term dominates in the high temperature 
limit and the combined compressibility approaches the 
non-interacting result as it should at very high temper- 
atures. At high temperatures the role of the exchange- 
correlation effects is diminished due to the increase of 
the thermal kinetic energy. As a consequence, the nega- 
tive compressibility at low densities becomes positive at 
a high enough temperature, T — Tp{'Kr'^/2y^^ . This be- 
havior (reversing sign due to increasing temperature) was 
observed in high mobility p-GaAs systemsi^. We have 
explicitly verified that our numerical HFA results agree 
precisely with the asymptotic high-temperature result of 
Eq. [T7]for t ^ 1. We note that the exchange correction 
to Kq/K falls off very slowly only as rg/ Vi ior t ^ 1, 
and as such quantum effects are quite large at large 
even for T 3> Tp. In fact, the classical regime in com- 
pressibility is approached only for T ^ r^Tp. 
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FIG. 1: (color online) (a) Calculated Ko/K in 2DEG as a 
function of rescaled temperature t for different values of rs. 
The inset presents Ko/K in the low temperature regime. The 
dashed and solid lines are corresponding to the analytical re- 
sults (given in Eq. [TS)) and the numerical results, respectively, 
(b) The inverse compressibility Kq /K in 2DEG as a function 
of Ts for different values of rescaled temperature t. (c) Calcu- 
lated tm, where the minimum of Kq/ K occurs, as a function 

of Ts. 



Numerical results of compressibility in 2D 
semiconductors 



In Fig.[TJa), we show the calculated inverse compress- 
ibility i^o / K as a. function of rescaled temperature t for 
six different values of r^. The red solid line corresponds 
to the non- interacting case {Kni) without the exchange 
interaction. As shown in Sec. III BI K^}{T] ~ h{t) t 



for i — )■ oo increases linearly in the high temperature 
classical limit. For finite the zero temperature inverse 
compressibility decreases with and eventually becomes 
negative (if > 7r/-\/2) due to the exchange energy. As 
shown in Fig.[lja) K^/K is not a monotonic function of 
rescaled temperature t. As t increases from zero K^^{t) 
initially decreases since the kinetic energy term is expo- 
nentially suppressed due to the Fermi surface restriction, 
but the exchange term keeps decreasing as t^\x\t. After 
reaching the minimum value at an rg-dependent charac- 
teristic temperature tm, increases linearly at high 
temperatures. This non-monotonic behavior becomes 
stronger at higher r^. More interestingly the negative 
value of at low temperature and high reverses 
its sign as the temperature increases. At high enough 
temperatures, where the kinetic energy dominates over 
the interaction energy, the system always has a positive 
compressibility, but typically strongly suppressed in mag- 
nitude from except for t — >■ oo. Thus, many-body 
effects manifest strongly in the compressibility even for 
T > Tp in most situations. 

In Fig. [ijb), we show our numerically calculated in- 
verse compressibility Kq/K as a function of rg for differ- 
ent values of the rescaled temperature t. It is clear that 
the interacting Kq/K manifests monotonically stronger 
many-body effects with increasing rg. The curves for 
t = and t = 0.8 cross each other at finite value of 
Tg, which corresponds to the non-monotonic tempera- 
ture dependence of inverse compressibility Kq/K in the 
low temperature regime. As an inset in Fig. [Ija) we 
explicitly show a quantitative comparison between our 
derived low-T asymptotic formula (Eq. I15p and the ex- 
act numerical HFA results verifying the presence of the 

In t term, which leads to the low-i minimum, in the HF 
Kq/K. In Fig.[IJc), we show as a function of our nu- 
merically calculated value of t{= T/Tp) — tm where the 
low-i minimum of Kq/K occurs. We note that although 
tm increases monotonically with increasing (decreasing) 
Ts (density), this increase is sublinear implying that the 
actual temperature (in Kelvin) Tm, where the minimum 
occurs, decreases (increases) with decreasing (increasing) 
carrier density n since Tp cx n and Vg oz n^^/^. This 
means that the experimental observation of this non- 
monotonicity of Kq / K as a function of temperature may 
be extremely difficult, if not impossible, with the experi- 
mental Kq / K manifesting only a monotonic increase with 
increasing temperature at all densities. The experimen- 
tal observation is further hampered by two additional 
complications: (1) the actual decrease in the magnitude 
of Kq/K associated with the shallow minimum is rather 
small (< 10%); (2) the low-density, large-rs regime of the 
2D system often develops strong disorder-driven density 
inhomogeneity. 

In Fig. [5] we compare our numerical results with our 
high-temperature analytic theory as given in Eq. [171 
The purely classical compressibility Kc{T) is given by 
Kq/Kc = t, which is the t{= T/Tp) oo limit of 
Kq/K,,, = h{t) = (1 - e-i/*)-\ as given in Eq. [IH We 
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FIG. 2: (color online) (a) Calculated Kq/K as a function 
of t for Ts = 1,4,8 comparing the HF result with the non- 
interacting result Ko/Kni ~ h{t) as in Eg. 1111 and the pure 
classical result Kq/Kc = t. Inset shows the HF result and the 
asymptotic form for the high temperature Kq/ K as in Eq. [17] 
with Ts = 8. (b) and (c) show K^i/ Khf for = 0.5, 1.0, 10, 
20 for t = ~ 20. 



note that the exchange correction to the non-interacting 
result, as given in Eq.[T7]falls off very slowly as r^i"^/^ ~ 
1/\/T. In Fig. Ilja) we show calculated Kq/K in the 
< = — 20 range for = 1,4,8 comparing the HF re- 
sult with the non-interacting result K^/Kni — h{t) as in 
Eq. [TTl and the pure classical result Kq/Kc = t. For the 
sake of comparison, we also show (for = 8) the corre- 
sponding high-temperature HFA result as in Eq. [T71 The 
interesting point to note here is that, as can be seen from 
Eq. I17[ the quantum exchange correction is quantita- 



tively substantial even at a temperature T ^ Tp for large 
rs with the HFA results being quantitatively well below 
(by - 20%) the classical result even for T = 20Tf. Al- 
though this appears somewhat counter-intuitive, the im- 
portance of quantum interaction persisting to high tem- 
peratures (T ^ Tp) for large can be understood by 
considering the relative magnitudes of the three dimen- 
sionless energy parameters = U/Ep, t = T/Tp^ and 
r = U/ksT that control the physics of an interacting 
quantum system - we note that rs, t, T are not indepen- 
dent parameters since t = rg/T. The classical nonintcr- 
acting limit requires both t, ^ 1 which necessitates 
t ^ as well. Thus one condition for the classical limit 
is i 3> when > 1 (and t 3> 1 for < 1), which is 
much stronger than i 3> 1 when is large! In Fig. 
and (c) we show Kni/KHP for r^. = 0.5, 1.0, 10, 20 for 
t = ~ 20 to emphasize that the exchange correction 
to the inverse compressibiUty is substantial in magnitude 
even for T > Tp. 



D. Results for 2D GaAs systems 

Since there has been considerable experimental 
activitjiii^ii^'i^ in measuring the 2D compressibility of 
both electrons^'^'^^ and holcs^^ in GaAs-based two- 
dimensional semiconductor systems, we provide in this 
section a set of numerical results for the HF compressibil- 
ity of 2D GaAs systems (electrons and holes) as a func- 
tion of carrier density (n in the unit of cm~^) and temper- 
ature (T in K). The parameters used in these numerical 
calculations are for electrons (holes): m = 0.07(0.4)TOe; 
g — 2 (i.e., = 2, g^ = 1); k = 13, where is the free 
electron mass in vacuum and k is the background lat- 
tice dielectric constant of GaAs-AlGaAs heterostructure. 
There is no valley degeneracy {gy = 1) in GaAs, and we 
consider the spin-degenerate {gs = 2) zero-magnetic field 
situation. 

In Fig. [S^a) and (b) we show our calculated dfj,/dn = 
{Kn^)~^ in the HFA for 2D GaAs electrons and holes 
respectively for T = 0,0.1,0.2,0.5,1,2,5,10,20,50,100 
K. The corresponding non-interacting T = result, 
d^i/dn EE [Do{Ef)]-^ = 0.24(TOe/TO)10-" meVcm^, is 
also shown as a constant horizontal line in each figure. 
The calculated temperature dependence is stronger for 
the holes than the electrons in Fig. [3] since the scale for 
the T-dependence, Tp, is much smaller (almost by a fac- 
tor of six for the same density) for the holes compared 
with the electrons because of the large difference in the 
two effective masses (0.07 versus 0.4). For n — 2 x 10^" 
cm-2: T| = 8.3 K, T^ = 1.5 K. Of course the qualitative 
behavior of K{n,T) is the same for both electrons and 
holes, it is only that the temperature scale for the holes 
is lower. 

The qualitative behavior of dfj,/dn as a function den- 
sity (n) and temperature (T) as shown in Fig. [3] is in 
very good agreement with experimental results^^, show- 
ing that the temperature dependence of compressibility 
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FIG. 3: (color online) (a) and (b) Calculated d^/dn = 
{Kv^)-^ in the HFA for 2D GaAs electrons (m = 
0.07me) and holes (m = QAme) respectively for T — 
0,0.1,0.2,0.5,1,2,5,10,20,50,100 K (from bottom to top). 
The solid horizontal line corresponds to the non-interacting 
T = result, dfi/dn = [Do{Ef)]~^. Insets of (a) and (b): 
Calculated rim, where the minimum of dji/dn occurs, as a 
function of T for electrons and holes respectively, (c) and (d) 
Zoom- in figures for the low-density parts of d^/dn versus n 
for 2D GaAs electrons and holes respectively. 



can indeed be very important in samples with low den- 
sities (and consequently with low Tp), particularly for 
holes because of their large effective mass. We note that 
Tp = 41.5 K (electrons) and 7.3 K (holes) for a density 
of n = 10^^ cm~^ (and Tp (x n). 

A particular qualitative feature of Fig. [3] deserves spe- 
cial attention: dfi/dn{(x K~^) as a function of n for a 
fixed temperature shows a very sharp upward turn with 
well-defined minimum at low densities, particularly for 



the low-T results. This striking low-n (and low-T) non- 
monotonicity arises from a quantum-classical crossover 
effect which turns on when Tp{n) < T from being 
Tp > T at some density and can therefore only be seen 
for low-T and low-n results. The minimum occurs at a 
density rim which we plot as a function of temperature 
T in the inset showing that Um increases with increasing 
T. To further emphasize this interesting behavior, we 
show in Fig.|3ljc) and (d), the low-density parts oi dfi/dn 
versus n for a few low T values where the minimum is 
clearly visible. This sharp increase of i^^^ or dfi/dn for 
low density (and low temperature) has been experimen- 
tally observed in 2D GaAs systems, and has often been 
associated with the 2D metal-insulator-transition (MIT) 
driven by disorder. By contrast, our theory does not in- 
clude disorder effects, we only include finite-temperature 
effects which are, of course, very strong at low densities 
where Tpin) < T. The intuition based on the T — 
theory of compressibility clearly must fail at some low 
density (i.e., for Tp{n) < T) since the T = limit of 
becomes large (in magnitude) and negative at very 
low n whereas the finite-T theory {for any finite T) pre- 
dicts that at the lowest densities, where Tp{n) < T, 
must become large and positive (i.e., the classical behav- 
ior for T ^ Tp), in fact eventually diverging as T/Tp 
for T ^ Tpl We believe, based on the results presented 
in Fig. [21 that experimentalists should re-investigate the 
older data for K^^ at low densities where finite temper- 
ature effects may be playing a significant role. There are 
some recent experimental results Ref. ^llj supporting our 
finding, but more experimental data are necessary in the 
interacting low density and low temperature regime to 
settle this question definitively. 

Given the possible qualitative importance of finite tem- 
perature corrections to the 2D compressibility as dis- 
cussed above (and shown in Fig. we provide in Fig. 2] 
the calculated Kq/K in the HFA for both electrons 
(Fig. HJa)) and holes (Fig. lUJb)) over a wide range of 
density and temperature. The non-monotonicity appar- 
ent for higher-T results in Figs. SJ^a) and (b) particularly 
for the hole data, is actually present in all the curves ex- 
cept that the non-monotonicity manifests itself at much 
lower density than the range covered in Fig. 2] for the 
curves that simply look like that Kq /K keeps on decreas- 
ing with density monotonically. In Figs. |4l[c) and (d) we 
show the high-T behavior (T = 20 - 100 K) of Kq/K. 

In Fig. [5l we show our calculated compressibility for 
2D GaAs electrons and holes as a function of tempera- 
ture for a few densities. Figs. [SJa) and (b) correspond 
to electrons showing Kq/K as a function of tempera- 
ture respectively over a wide T = — 100 K range [IJa)] 
and a narrow low temperature range T = — 5 K ^h)] 
whereas Figs. [5]Jc) and (d) show the same for 2D holes 
in GaAs. In Figs. Efb) and (d), we show a comparison 
between our low-T analytical and numerical results also. 
The important point to note in Fig. [5] is that the low- 
temperature minimum should be observable in 2D semi- 
conductor systems in careful measurements of Kq/K as 
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FIG. 4: (color online) (a) and (b) Calculated Ko/K as a 
function of carrier density in the HFA for electrons (m = 
0.07me) and holes (m = 0.4me), respectively, over a wide 
range of temperature T = — 10 K. (c) and (d) show the 
high-T behavior (T = 20 - 100 K) of Kq/K for electrons and 
holes respectively. The temperature increases from bottom to 
top in each figure. 



a function of temperature provided disorder effects are 
unimportant, i.e., highest- mobility samples are used for 
the study. Lower density samples would typically mani- 
fest deeper minima as can be seen in Fig. [SJ 



III. COMPRESSIBILITY IN GRAPHENE 

In this section, we theoretically calculate the finite 
temperature compressibility of monolayer graphene in- 
cluding exchange interaction effects. It has been shown 
that the HFA is an excellent quantitative approximation 
due to the small contribution of the correlation energy in 
graphene^''. We focus on extrinsic graphene, i.e., gated 
or doped graphene with a tunable 2D free carrier density 
n of electrons (holes) in the conduction (valence) band, 
i.e., the chemical potential /i being positive (negative). 
For undoped (intrinsic) graphene the chemical potential 
is zero even at finite temperatures, which gives rise to 
a logarithmically divergent compressibility if disorder ef- 
fects are neglectecU^iii. The basic feature of graphene 
compressibility and its functional dependence on density 
have been well studied in the literature at T = 0^"^'-^. 

We use the same definition for the compressibility 
given in Sec. |TT1 However, there are two main differences 
between the usual 2D semiconductor and graphene. The 
first one is the difference between their energy dispersion 
relation: 2D semiconductor has parabolic dispersion re- 
lation while graphene has linear dispersion. The kinetic 
energy in graphene is given by ej^^s = svF\k\, where |fc| 
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T (K) T (K) 

FIG. 5: (color online) Calculated Ko/K in the HFA as a 
function of temperature for various densities. Figs, (a) and (b) 
are for electrons(m = 0.07me), (c) and (d) are for holes(m = 
0.4me). (a) and (c) show Kq/K in the HFA for a wide range 
of temperature (T = — 100 K) for various densities n = 
2, 4, 6, 8, 10, 12, 14 X 10^° cm"^ (from top to bottom), 
(b) and (d) show Kq/K at low temperatures (T = — 5 
K) for various densities n = 2, 4, 6, 8, 10, 12, 14 x 10^° 
cm~^ (from bottom to top) and compare the low-T analytical 
results (dashed lines) and numerical results (solid lines). 



is the wave vector with respect to the Dirac point, vp 
is the Fermi velocity with the value oi vp ^ 10^ cm/s 
and s = ±1 for the conduction and valence band respec- 
tively. The second difference is that we need to consider 
the contribution from valence band electrons in graphene 
because graphene is a gapless semiconductor so that the 
valence band and the conduction band touch each other 
at the Dirac point. In graphene the total degeneracy 
(7 = 4 (from spin degeneracy gs — '2 and valley degener- 
acy ffi, = 2). 

Since the non-interacting chemical potential of 
graphene is fj,o = Ep = hvpkp at zero temperature the 
corresponding non-interacting inverse compressibility is 
given by 

^n'^TThvp/{2kp), (18) 

where we use g = ^ and the electron density is related to 
the Fermi wave vector as = irn. Unlike the regular 
2D systems, noninteracting d/x/dn in graphene is density 
dependent (i.e. d^i/ dn oc n~^/^) due to the linear energy 
dispersion of graphene. 

The chemical potential within HFA is the sum of the 
non-interacting kinetic energy part and the exchange self- 
energy part. The exchange self-energy for graphene is 



9 



given a! 



13. 



E 



—-V{\k-k'\)nF 
(27r)^ 

(19) 

where s,s' = ±1 indicates the band indices, V (q) — 
/ {nq) is the bare Coulomb interaction (k is the 
background dielectric constant), and Fss'{k,k') = (1 + 
ss' cos 9 kk>)/ 2 arises from the wavefunction overlap fac- 
tor where Okk' is the angle between fe and k' . In regular 
(i.e., non-chiral) 2D systems, F = 1, but in graphene, 
which is a chiral material, the chiral factor F is impor- 
tant due to the underlying pseudospin dynamics. With 
Eq. we calculate the interacting compressibility of 
graphene at finite temperatures. Since the zero tempera- 
ture exchange self energy can be calculated analytically^'^ 
we first consider the zero temperature compressibility as 
the starting point of our discussion. 



A. Zero temperature compressibility 

At r = 0, the Fermi distribution function in Eq. [TOl 
becomes npiC} = Q{Ef — 0- We divide the exchange 
self-energy into two parts^^: one from the intraband tran- 
sition E!jf', and the other from the interband transition 



That is, Ex,3(fc) = T.'^Mk) 



T.f:l{k), where 



dq 



(2^) 



\k,k~ q); 



Vcxt 



(fc) 



/ 7T7^'5"F(Cfc-q,s')K(9) F,s'{k, k-q) 

(20) 

where 5nF{^k-q,s') = npi^k-q^s')- ^{'^- s') is the differ- 
ence in the electron occupation from the intrinsic T = 
case. 

After some algebra (the detail of the derivation is given 
in Ref. [11]), we find the density dependent exchange 
contribution to the inverse compressibility Kq /K at zero 
temperature 



K 



1 



iln^4 

4 \y 







1 

2^ 



(21) 



where y = kp/kc with a momentum cut-off /cc ~ a ^ 
where a is a lattice constant of graphene, C w 0.916 is 

2 

Catalan's constant, and r. 



is the graphene cou- 



pling constant (or fine structure constant). We note that 
the Iny = lii(kp/kc) term takes care of the divergent 
compressibility of intrinsic graphene which does not en- 
ter our discussion in any significant manner. Unlike in 2D 
systems, where the interaction parameter Ts ^ rt^^/^, for 
graphene is a constant in density due to its linear en- 
ergy dispersion. However, by adjusting the background 
dielectric constant (k) we can vary the value, from 
Ts = 2.2 (for K = 1, graphene suspended in vacuum) 



to very small by making k very large. For graphene 
on Si02, rs = 0.8. Note that unlike ordinary 2D sys- 
tems the calculated graphene compressibility with the 
exchange correction is always positive, which has recently 
been measured by several different techniques recentlj*ii^. 



B. Finite temperature compressibility 

In this subsection, we present the theoretical formal- 
ism of the finite temperature Kq/K and its asymptotic 
analytical formula at low temperatures (t <C 1). We con- 
sider extrinsic graphene with Ep > and concentrate 
on the situation with the chemical potential lying in the 
conduction band with no loss of generality. 

The finite temperature chemical potential (without ex- 
change energy) /j,o must be calculated by the conservation 
of the total electron density, i.e. n(T) = n{T = 0)+p{T), 
where n(T) and p{T) are the electron and hole density 
at T, respectively. They are given by 



f d^P 



1 



(27r)2 exp[/3(£p - m)] + 1 



and 



p{T) = g 



d^p 



(27r)2 exp[-/3(£p - /i)] + 1 ' 



(22) 



(23) 



where g is the total degeneracy and f3 = l/ZcsT. Thus, 
we have the self consistent equation for fi as 



1 /Tp 



2 \ T 



- Fi iPfi) - Fi i-Pfi) , (24) 



where Fn{x) is given by 

Fn{x) - 



V'dt 



1 -|- exp(i — x) 



(25) 



Then we obtain the non-interacting chemical potential 
for both low and high temperature limits for graphene as 



A^o(T) 
Mo(T) 



Ep 



2 / rj-i 



Ep Tp 
41n2^ 



for T/Tp < 1 (26) 
for T/Tp > 1.(27) 



With Eqs. ([26| and (p7)) we can find the noninteract- 
ing graphene compressibility at low temperatures and at 
high temperatures, respectively. Let no/Ep = /Iq and 
t = T/Tp^ then the non-interacting part of the inverse 
compressibility can be written as 



Kg 
Km{T) 



TT t 
1 1 



21n2 t 



fort<l, (28) 
for i > 1. (29) 
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Thus the leading order correction to the inverse com- 
pressibihty Kni{T)~^ increases quadratically at low tem- 
peratures and then it decreases inverse linearly at high 
temperatures. This fact indicates that the graphene com- 
pressibility shows the opposite behavior to the ordinary 
2D systems as shown in Sec. HIl where the inverse com- 
pressibility of 2D systems decreases first and then in- 
creases as the temperature increases. We emphasize that 
this qualitative difference in the noninteracting finite- 
T compressibility between 2D systems and monolayer 
graphene arises both from the linear, chiral energy dis- 
persion of graphene and from its gaplessness allowing for 
thermal excitations of conduction valence band electron- 
hole pairs at finite T. 

Now we calculate the full compressibility at low tem- 
peratures, including the exchange energy. As done for 
the zero temperature case, the exchange self energy of 
graphene can be separated into contributions from the 
intrinsic part S"'' and the extrinsic part S™'^^. We can 
see that the first term in Eq. (15(1)) is E"'* while the second 
term is For /i > the Hartree-Fock self energy in 

graphene at finite temperature becomes 



1 — COS0 



E^"'(fc,r) 27rK -^-^ o'^P('3(-'ii'-F9-M)) + l 2 



(30) 



qdqdO 



l-Hcosg 



E°''*(fc,r) 27rK \k-q\ cxp(^;(toFg-po))+i 2 

The intrinsic exchange self-energy E!,,"' in Eq. (pO)) can 
be derived by using the fact that the exponential term in 
E"* is exponentially suppressed and almost equal to 1 in 
the low temperature limit. Then, we have^^ 



E-*(r) 



E{y) 



TT , ,4, 71 

h - In(-) 

E{x) - -^x") 



0(e 



(31) 

where y = kp/kc ^ 1, Ki{x) and E{x) are the complete 
elliptic integral of the first and second kinds, respectively. 
Note that the leading order temperature correction of the 
intrinsic exchange self-energy at low temperatures is ex- 
ponentially suppressed. On the other hand, the extrinsic 
part of the exchange self-energy becomes 



Er(T) 



e^kp 



C 



1 



TT^t^ Int 



1_^, 

2 12 ' " 12 

(32) 

0.916 is Catalan's constant, = — 1 + 27 — 
s 0.877 where G ~ 1.282 is Glaisher's 
0.577 is Euler's constant. Differen- 



where C ? 
121n(G)+ln(if) 
constant and 7 :i 
tiating the total chemical potential with respect to the 
density, we have the asymptotic form of inverse compress- 
ibility in the low temperature limit 



Kn 



K{T) 



= fit) + 



4 Vw/ 8 TT 2n 



12 



(33) 
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FIG. 6: (color online), (a) Calculated Kq/K in graphene as a 
function of temperature t for three different values of density 
in the units of lO^^cm"^. (b) Calculated Ko/K in graphene 
as a function of carrier density n for six different values of 
rescaled temperature t. In this calculation the momentum 
cut-off kc = l/a {a = 2.46 A) and k — 2.5 are used, (c) The 
semi-log plot of Kq/K versus carrier density with different 
temperatures T for n = 1 and (d) for k = 2.5. 



where f{t) = 1-1- is the non-interacting part of the 
inverse compressibility given in Eq. ()28|) . In the low tem- 
perature limit the leading order temperature correction 
to the total inverse compressibility is the same as that 
in the 2D semiconductor system due to the exchange en- 
ergy, i.e., t^lnt. However, the logarithmic correction can 
be detected only at very low temperatures {t < 0.01). In 
general, the correction from the non-interacting part 
dominates at low temperatures {t > 0.01) and the posi- 
tive coefficient of t'^ term gives rise to increasing behav- 
ior of K^^ with temperature. Our asymptotic results at 
low temperatures agree very well with the full numer- 
ical results shown in Fig. [6l[a) . The results provided in 
this section generalize the existing graphene literature on 
compressibilityi^ to finite temperatures. 



C. Numerical results of graphene compressibility 

In this subsection, we present our numerical results of 
finite temperature compressibility in graphene. Through- 
out this section, we use vp = 10* cm/s and the wave 
vector cut-off kc = 1/a {a — 2.46 A). We show re- 
sults for suspended graphene (k ~ 1.0) and for graphene 
on Si02 {k ^ 2.5). In Fig. [6ja) we show the tem- 
perature dependence of the total inverse compressibil- 
ity Kq/K for three different values of carrier density 
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FIG. 7: (color online), (a) Calculated Ko/K in monolayer 
graphene as a function of temperature t for three different 
values of background dielectric constant within HFA for car- 
rier density n — lO^'^cm"^. The dot-dashed (orange) line 
presents the temperature dependence of the non-interacting 
Ko/K. (b) Calculated Ko/K in monolayer graphene as a 
function of the background dielectric constant k for a few 
fixed values of t. kc — 1/a {a — 2.46A). 



cm^^ with dielectric constant k — 2.5 (corresponding to 
graphene on the Si/Si02 substrate). We can clearly see 
that the temperature dependence is non-monotonic and 
this non-monotonic behavior mostly comes from the non- 
interacting compressibility, which is different from the 
2D semiconductor case. Unlike the 2D parabolic-band 
case, the compressibility in graphene does not change 
sign in the range of experimentally relevant parameters 
and typically exchange corrections are always quantita- 
tively small (< 20%). We show the carrier density depen- 
dence of Kq/K in Fig.[6jb), where the Kq/K is a mono- 
tonically decreasing function of carrier density for small 
values of t (T ~ 400K for i = 1 and n = 0.1 x IQi^cm-^). 
In Fig. |ni[c) and (d), we present the calculated Kq/K as 
a function of carrier density for different temperatures. 
For fixed temperatures, Kq/K is a decreasing function of 
n at higher carrier density, which corresponds to the low 
temperature limit manifesting n^^ density dependence. 
On the other hand, at lower carrier density, Kq/K has 
n^^^ density dependence, which corresponds to the high 
temperature limit and the main contribution to Kq/K 
comes from the non- interacting part (see Eg. [^^ . 

In Fig. [71 we compare the temperature dependence for 
different values of background dielectric constant k ~ 1 
(i.e. graphene in vacuum), k = 2.5 (graphene on Si/Si02) 
and K = 15 (graphene on hafnium oxide Hf02), respec- 
tively. The value of represents the fine structure con- 
stant of graphene which depends only on k and does not 
depend on the carrier density because of the linear en- 
ergy dispersion. Larger dielectric constant corresponds 
to smaller values of Ts, which indicates a weak-coupling 
system in terms of electron-electron interaction^^. Note 
that the trends of Kq/K among different values of di- 
electric constant are similar. The larger the value of di- 
electric constant (the smaller the value of r^), the less 
strong is the density dependence of Kq/K. We also 
show in this figure our calculated Kq/K in HFA com- 



pared with the noninteracting Kq/K as a function of 
the background dielectric constant k for a few fixed val- 
ues of t. We could see that the inverse compressibil- 
ity Kq/K first increases as the temperature t increases, 
reaches a maximum, and then decreases as the tempera- 
ture t further goes up. This non-monotonic behavior of 
Kq/K dependence on the temperature is mainly due to 
the non-interacting part of the chemical potential, which 
is consistent with the calculation of the non-interacting 
part of the chemical potential (Eal?7| . From Eq (l?fl) we 
know that the temperature dependence of the chemical 
potential is (with a positive slope) at low tempera- 
tures while the temperature dependence of the chemi- 
cal potential approximates as 1/t (with a negative slope) 
at higher temperatures. There must be an intermediate 
point at which the chemical potential changes its sign of 
slope and so does the compressibility K . Since K^^ is 
the derivative of chemical potential with respect to n, it 
has similar temperature dependence as the chemical po- 
tential. Therefore, the non-interacting part of Kq / K be- 
haves as at low temperatures and then changes to \/t 
behavior. The domination of graphene compressibility 
by noninteracting effects is just a direct manifestation of 
the weakness of electron-electron interaction in graphene 
compared with the 2D semiconductor case. 



IV. COMPRESSIBILITY OF BLG 

In this section, we calculate the BLG compressibil- 
ity within the HF approximation. There is some re- 
lated theoretical work in the literature dealing with the 
problem of compressibility in 'SLG^^^^^— but the spe- 
cific and detailed temperature dependent compressibil- 
ity presented in the current work is not available in 
the literature. To calculate BLG compressibility at fi- 
nite temperatures we use the two-band approximation, 
which is valid for the low density limit. At high den- 
sity, BLG dispersion approaches the linear energy dis- 
persion of monolayer graphene as the Fermi energy be- 
comes large^^, and therefore the high-density BLG com- 
pressibility behaves similar to the MLG compressibility 
studied in the last section. It has been shown that the 
two band model presents the most important qualita- 
tive signatures of the BLG compressibilityJ^iiS'^Sii^ For 
low density bilayer graphene within two-band model, the 
quadratic approximation is commonly used, i.e., ^fc = 
h'^Vpk'^/'ji = h^k^ /{2m), where 71 is the interlayer hop- 
ping parameter from the tight-binding approximation 
and m = ^i/{2v\) is the effective mass of the electrons. 
This dispersion comes from a low-energy effective the- 
ory of bilayer graphene^, which essentially discards the 
two split bands and confines electrons to those lattice 
sites not involved in the interlayer coupling. While the 
quadratic dispersion is relevant at low energies, the ac- 
tual dispersion of the bands is hyperbolic. At large wave 
vectors, relevant at large densities, the dispersion is ef- 
fectively linear and the system should behave like single 
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layer graphene. Thus our results obtained in the previous 
section should be applicable to the high density BLG. 
The exchange self energy of BLG is given by 

Sx.s(fc) = -5m^k-k'riF(a',.')^.v(k,k'), (34) 



where Vk_k' — 27re^/K|k— k'| is the Coulomb potential, 
= fc^/2m, npix) is the Fermi distribution function, 
and Fs_5/(k, k') = (1 + ss' cos(26'kk'))/2 is the wave func- 
tion overlap factor and 9kk' is the angle between k and k'. 
To get the explicit temperature dependence we rewrite 
Eq. dMl) as 



+ i^Vk_k'[/(^k') + .9(a')][l -cos(20)p5) 

k' 

where f{x) = l/(e'^(^-'^) + 1) and gix) = l/(e''(^+^) + 1), 
^2D{k) is the ordinary 2D self energy given in Eq. (5), 
and S™* is given by 



d^k' 1 - cos(2g) 
^14-. (36) 



3ttk 



-~p{z) + q{z) 



where 



p{z) 



and 



3z3 

(l + 7z^)g(l/z) 
3z2 

(l + 2z2-3z^)j^i(z) 

3^3 

4(z 2-l)Xi(l/z) 
3z2 



z < 1 



z > 1 



(37) 



z < 1 



z > 1 



(38) 



Here z = k/kc, kc ^ 1/ a is the momentum cut-off (sim- 
ilar to the MLG situation) with a lattice constant a. 
Ki{x) and E{x) are the complete elliptic integral of the 
first and the second kind, respectively. Since we have 
studied S2D in Sec. II and S™* is independent of tem- 
perature, the new feature of temperature dependent com- 
pressibility of BLG arises from the third term in Eq. ([35]). 



^i'Hfc) = lj2^k-Af{^k')+giCk'm-cos{2e)]. (39) 



Then the total temperature dependent HFA chemical po- 
tential of BLG can be calculated to be 



(40) 



A. Zero temperature BLG compressibility 

At T 0, using f{x) = 61(1 - x) and 5(2;) = in 
Eq. ([M)) . the exchange self-energy can be calculated as: 



Sx(fc) 
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TTK V kp 
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2e^k r 



37rK 
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q 



Then, the exchange self-energy at A: = fc^ is given by 



kp , 
(41) 



'^x{kF 



2e^kF 



2e^kF 



p{y) + qiy) . (42) 



where y = kp/kc- The total chemical potential for BLG 
within HFA is then written as 



Ei. 



1-^ (^ + ^b(y)-9(y)] 



(43) 



where = 2e'^m/ {nkp) is the dimensionless BLG cou- 
pling. In the following numerical calculation, we use the 
the value of dielectric constant k = 2.5 and effective mass 
m = 0.03toji^. Differentiating Eq. |321 with respect to n 
and making use of the relation k"^ = tto, we get the zero 
temperature BLG inverse compressibility within HFA 

^ = 1 - 7 + \[p{y) - q{y)]) + ^^(y) (44) 

where is the no n- interacting compressibility of BLG, 
given by 



and 



i^ + y 



1 -KIT? 



2/2 _i 



(45) 



w{y) = < 



3y3 ^(y) + ^^^^^(y) 



l + y2 



1-r 



1^ E{l/y) + ^^K^{l/y) y>l 

^ ^ (46) 

At low enough densities the zero temperature compress- 
ibility within HE approximation becomes negative simi- 
lar to the corresponding 2D semiconductor system. The 
corresponding density is ric ^ 10^" cm~'^ or Vs ^ 27, for 
which the negative compressibility manifests itself. This 
is a very small (large) value of density (r^) compared 
with the ordinary 2D case where the negative compress- 
ibility shows up for Ts > 2.2. However the negative 
compressibility obtained in HFA has not been experi- 
mentally observed in bilayer graphene systemsS^"— . The 
non-negative compressibility may be attributed to cor- 
relation effects. It is shown that the contributions from 
exchange and correlation almost cancel each other out in 
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FIG. 8: (color online) (a) Calculated Kq/K in BLG as a func- 
tion of temperature t — T/Tf for different values of r^. The 
inset presents Kq/K for = 2.0 showing the non-monotonic 
behavior, kc = 1/a {a = 2.46 A) and k — 2.5 are used, (b) 
The inverse compressibility Kq/K in BLG as a function of 
Vs for different values of temperature t. (c) d^/dn in BLG 
versus carrier density n for various temperatures. Inset shows 
the details in the low density limit in the semi-log scale. 



B. Finite temperature BLG compressibility 

At finite temperatures the non-interacting chemical 
potential is temperature independent in BLG, i.e., 



^1Q{T) = EF. 



(47) 



Then the non-interacting compressibility of BLG is also 
temperature independent and is the same as the zero 
temperature compressibility given in Eq. (|45p 



(48) 



At finite temperatures, the temperature dependence of 
BLG compressibility comes from the exchange contribu- 
tion. This is very intriguing because the non-interacting 
kinetic energy dominates exchange self energy at high 
temperatures for both ordinary 2D systems and MLG. 
But in BLG the entire temperature dependence arises 
from the exchange energy! Thus, any experimentally ob- 
served temperature dependence in the BLG compressibil- 
ity must entirely be a many-body effect, at least within 
the 2-band approximation (i.e., at not-too-high carrier 
densities). 

Let us consider the exchange self energy, Eq. (I55t . at 
finite temperatures. The first term in Eq. ([35]) has been 
calculated in Sec. [ll] [see Eq. (|14l) ] and the second term 
is temperature independent. The third term [Eq. (j39l) ] 
becomes in the low temperature limit t = T /Tp <^ 1: 



Sttk 



.(£ _^_!Lt2) + 0(t4). (49) 



Thus, we have the total HFA chemical potential in the 
low temperature limit 



m(t) 



1 



TT 



TT 

Be 
4 



2j-2 



72 



1 2 2, 16 
48 t 



(50) 



where B = [ln(47r) + 12C'(-1)] ~ 0.449 and C'(-l) is 
the derivative of the zeta function. Differentiating the 
asymptotic formula of chemical potential with respect to 
the density, we have the BLG Kq/K for i <C 1 



K{T) ~ K{T : 



0) TT 



0.55i2 + —tHnt 



(51) 



BLG at low density, leaving a small positive compressibil- 
ity at all densities in contrast to the corresponding 2DEG 
systemsir— where the electronic compressibility can be- 
come negative at low densities.— 1^ The non-quadratic 
band structure also gives rise to significant affects on 
the compressibility as has recently been discussed in the 
literature^ii 



where Kq/K{T — 0) has been derived in Eq. We 
have the same leading order term In t arising from the 
exchange self-energy as the ordinary 2D system. Thus, at 
low temperatures the calculated inverse compressibility 
decreases as the temperature increases. 

In the high temperature limit T 3> If, the chemical 
potential can be calculated to be: 



ti{T) = Ep 



(52) 
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where Bi = [-Sv^tt^ + 8 (-1 + 2^2) 0FC(|)]/64 « 0.4 
and Cix) is the zeta function. Then the inverse compress- 
ibihty is given by 



(53) 



In the high temperature hmit (t ^ 1) the inverse com- 
pressibihty of the ordinary 2D systems increases hnearly 
with temperature and approaches the non-interacting 
value. However, since the non-interacting compressibihty 
of BLG is temperature independent and the most dom- 
inant exchange term decreases as the high tem- 
perature inverse compressibihty of BLG approaches the 
following high-temperature limit: 



Ko 



K{T) 



1 



1 a 

8 as ' 



(54) 



where = T? k j {me^') is the BLG effective Bohr radius. 
Since a jag < 1 the negative inverse compressibility re- 
verses its sign to the positive values at high temperatures 
and asymptotically approaches a smaller value than the 
zero temperature inverse compressibility. This intrigu- 
ing result is purely a high-temperature manifestation of 
exchange effect within the 2-band BLG approximation. 

In Fig. ISl we present our numerical results of tem- 
perature dependent inverse compressibility in BLG sys- 
tems. In the calculation we use = 1/a (a = 2.46A) 
and K = 2.5. As shown in Fig. [8] the calculated 
BLG inverse compressibility shows very weak tempera- 
ture dependence. The weak non-monotonic behavior of 
temperature-dependent inverse compressibility is shown 
in the inset of Fig. ^a). In Fig. El^b), we present the 
inverse compressibility versus r^. decreases mono- 

tonically with Ts and becomes negative at 27. In 

Fig. I3c)) we show dfi/dn as a function of density and 
the inset shows the same figure in the semi-logarithm 
scale to show clearly the behavior at low densities. At 
densities n < 10^'^cm~'^ the zero temperature dfi/dn is 
negative, but as temperature increases it reverses its sign 
at < sa [0.4/(7r/rs — o/Sas)]^. As temperature increases 
further it approaches (fi.^7r/2m)(l — a/Sas) regardless of 
density. 

In Fig. [9l we show, for the purpose of comparison, the 
calculated Kq/K for the 2D, MLG, and BLG systems as 
functions of and t in order to demonstrate the quali- 
tative difference between the results. 



V. DISCUSSION AND CONCLUSION 

In this work, we have calculated the finite tempera- 
ture compressibility of graphene and ordinary 2D semi- 
conductor systems within the Hartree-Fock approxima- 
tion. We present both analytical and numerical results 
of Kq/K as a function of temperature, the dimension- 
less interaction parameter r^, and density. We find that 
the calculated temperature dependent compressibility is 
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FIG. 9: (color online) Calculated Kq/K for the 2D, MLG, 
and BLG systems as a function of t for a few values of r^. 



non-monotonic in both graphene and 2D semiconductor 
systems. In monolayer graphene, the inverse compress- 
ibility K~^ increases with temperature at low temper- 
atures and decreases at high temperatures, reaching a 
maximum value at an intermediate temperature. This 
non-monotonicity arises entirely from the behavior of the 
non-interacting graphene compressibility. The tempera- 
ture dependent inverse compressibility in 2D semiconduc- 
tor systems decreases with temperature at low tempera- 
tures and increases at high temperatures. In BLG the in- 
verse compressibility decreases as temperature increases 
at low temperatures, but it approaches a value which is 
less than the non-interacting value at high temperatures. 
The leading order temperature correction to the inverse 
compressibility K~^ in HFA for both graphene and 2D 
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semiconductor systems is Int. Our analytic results are 
in agreement with our numerical results both at low and 
at high temperatures. 

Our use of the Hartree-Fock approximation is not a 
particularly restrictive approximation for the theory of 
the electronic compressibility since it is well-known that 
Hartree-Fock theory works well for the calculation of the 
compressibility and its density dependence, at least at 
T = 0. The Hartree-Fock approximation gives results 
numerically very close to the full dynamical RPA for the 
2D compressibility at T = 0^^, and the Hartree-Fock 
compressibility results are in good agreement with the 
compressibility measurements for 2D semiconductors^ 
and graphene^. In fact, the Hartree-Fock approxi- 
mation for 2D semiconductors is in better quantitative 
agreement with the density dependence of the measured 
compressibility^ than the corresponding RPA theory as 
long as the strict 2D approximation is used since the fi- 
nite width corrections for the realistic 2D semiconductors 
tend to cancel out the correlation corrections neglected 
in the HFA. In any case, we consider the HFA as the first 
step necessary for understanding the temperature depen- 
dent 2D compressibility in interacting electron systems, 
and only future work, particularly experiments compar- 
ing with our predictions, can establish the necessity of 
improved theoretical treatments involving dynamical cor- 
relations neglected in the HFA. 

Our main finding that the temperature dependence of 
compressibility is much stronger in 2D systems than in 
graphene is consistent with existing experimental results. 
In particular, our results of Sec.|ll]agree well with the ex- 
perimental study of Shapira et al.^^ on 2D GaAs holes 
who discovered substantial temperature dependence in 
the 2D hole compressibility in the 1 — 100 K temperature 
regime with the experimental data showing reasonable 
qualitative agreement with our theoretical HFA results 
presented in Sec. [Tib. Detailed quantitative comparison 
with the experimental data of Shapira et al.,^^ is not 
particularly meaningful (and is not attempted) because 
of a number of reasons including our use of the exchange- 
only Hartree-Fock approximation, our neglect of the fi- 
nite thickness of the 2D semiconductor structures which 
is often important for quantitative considerations, our ne- 
glect of disorder effects which are certainly important at 
lower carrier densities, and the lack of information about 
some essential experimental parameters (e.g., depletion 
charge density). 

In addition to the experimental study of the tem- 
perature dependent 2D hole compressibility by Shapira 
et al.,^^ discussed above, there have been several low- 
density studies of the temperature dependence of the 
compressibility in both electron and hole 2D systems-'^'iS 
carried out in the context of 2D MIT (i.e., disorder- 
induced density-tuned localization of 2D semiconduc- 
tor systems). As we discussed in Sec. [Hi these low- 
density experimental investigations should be revisited 
in light of our current work demonstrating the impor- 
tance of the temperature dependence of Kq/K in the 



low carrier density regime by virtue of the dimension- 
less temperature T/Tp being large precisely in the low- 
density regime (where disorder effects are also strong) 
since Tp cx n in 2D systems. Unfortunately, tempera- 
ture, disorder, and interaction effects are all strong in 
the low density regime, making any quantitatively reli- 
able theoretical work essentially impossible in the low 
carrier density regime, and our work establishes the im- 
portance of finite temperature effects on the compress- 
ibility in the low carrier density regime. In particular, 
the density inhomogeneity and puddle formation that 
happens at low carrier densitj*21i^ due to the failure of 
screening would very much complicate the observation 
of the low-density (and low-temperature) features in the 
compressibility predicted in our theory unless one uses 
extremely high-quality samples with very low disorder. 
Recently, there has been a theoretical investigation^'' of 
the effect of disorder on low-density graphene compress- 
ibility clearly establishing the importance of disorder in 
low-density compressibility. 

Since the Fermi temperature Tp goes as Tp(iii K)= 
4.2n (2D GaAs electrons); 0.74n (2D GaAs holes); 136v^ 
(monolayer graphene); 4.2n (bilayer graphene), where 
n = is the carrier density measured in the units 

of 10^° cm~^, it is obvious that the quantitative effect of 
finite temperature, even at low carrier densities, is by 
far the strongest (weakest) in 2D GaAs holes (mono- 
layer graphene), which is consistent with experimental 
observations. Since disorder effects become important for 
n Ui, where rii is the background random charged im- 
purity density in the environment, the low density regime 
associated with large T/Tp values becomes even more 
challenging to achieve in the laboratory experiments on 
graphene. We expect our predicted temperature depen- 
dence to manifest in graphene compressibility in room 
temperature experiments in very high-mobility samples 
where T/Tp ^ 1 values may be achievable. 

An important relevant question for our theory is: what 
is the most suitable system for the experimental observa- 
tion of our theoretical predictions? Obviously, our pre- 
dicted high-temperature behavior, where the theory is in 
very firm ground since the exchange energy is likely to 
be the exact leading-order many-body correction to the 
non-interacting compressibility in the high-temperature 
limit, should be valid for all systems and should be 
observable in 2D semiconductors (both electrons and 
holes) rather easily by measuring the compressibility for 
T = 10- 200 K for electron (hole) densities- 10^° - 10" 
cm~^ so that T > r^Tp condition is satisfied. Given 
that Tp < Tp at the same density, the high-temperature 
behavior is much more easily observable in 2D GaAs 
holes, and has, in fact, already been observed at least 
in one experimental study in 2D holes^'^. While 2D 
GaAs holes are the obvious ideal candidates for observ- 
ing our predicted high-temperature compressibility be- 
havior, monolayer graphene may turn out to be not par- 
ticularly well-suited for the temperature-dependent com- 
pressibility studies because of its very high relative Fermi- 
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temperature {Tp'^'^ ^ 1350 K for a doping density of 
10^^ cm"^) and its relatively weak temperature depen- 
dence. Obviously, the condition for observing the high- 
temperature behavior of compressibihty is the abiUty to 
reach T > Tp which necessitates lower Fermi tempera- 
tures. By contrast, the low-T behavior associated with 
the In T term in the compressibility may be better ob- 
served in a system with a relatively high value of Tp so 
that the regime T/Tp < 1 can be explored over a fairly 
broad range of temperature with the temperature tm not 
being too low. Thus, very high mobility n-GaAs het- 
erostructures or quantum wells with n ^ 2 x 10^° cm~^ 
so that Tp ~ 8.3 K (and Tm ~ 3 K) may be the ideal 
system to search for our predicted shallow minimum in 
Kq/K. One serious problem is that the minimum may 
be too shallow to be uniquely determined experimentally. 

Finally, we discuss our very interesting analytical find- 
ing of the non-monotonicity in graphene inverse com- 
pressibility associated with the T^ In T correction we ob- 
tain analytically (and verify numerically). Although this 
would not be an easy effect to detect experimentally be- 
cause of its quantitative weakness (and became the sub- 
leading correction goes as and thus the ratio of the 
two, ln(r/rF'), is always challenging to observe in ex- 
periments even under the best of circumstances), it is 
nevertheless interesting to discuss its origin and its ro- 
bustness beyond our approximation scheme. Our results 
are exact within the Hartree-Fock approximation scheme, 
and there is no doubt that the exchange self-energy con- 
tributes a In r leading-order contribution to the in- 
verse compressibility. The first question is whether this 
is intrinsically a dimensionality effect occurring only in 
two dimensions. We have therefore carried out the corre- 
sponding analytical calculations for the 3D temperature- 
dependent compressibility (presented in the Appendix |A] 
of this paper), finding that Kq/K has a T^lnT contribu- 
tion in three dimensions also. Thus, it appears that the 
interesting non-monotonicity associated with the In T 
correction is a results of the Hartree-Fock approximation 
arising from the exchange self-energy diagram, and is not 
intrinsically a two-dimensional effect. 

We have investigated this question by calculating the 
interacting compressibility in the screened Hartree-Fock 
approximation (sometimes also called static RPA) for the 
2D system, where the bare Coulomb interaction 'F' ap- 
pearing in the exchange self-energy is screened by the 
static RPA dielectric function. (These results are shown 
in Appendix [Bl) The screened HFA (Appendix [B]) or 
static RPA does not contain the In T term in the finite 
temperature inverse compressibility, but is not a reliable 
approximation at all since at T = it predicts that Kq j K 
should be always positive for all values of in clear dis- 
agreement with experimental finding^'^. Thus the ab- 
sence of the In T term in the screened HF approxima- 
tion cannot be taken seriously since the corresponding 
T = result is in qualitative disagreement with experi- 
mental results. It is interesting to speculate whether the 
T^ In T term survives higher-order diagrams associated 



with dynamical correlation effects, and this remains a 
challenge for the future. Although it is fairly straight- 
forward to calculate the compressibility including corre- 
lation effects at T = — for example, the ring-diagram 
contributions to the compressibility can be exactly calcu- 
lated at T = — it is a formidable challenge, both nu- 
merically and analytically, to extend to the correspond- 
ing ring diagram calculations of interacting compressibil- 
ity to finite temperatures. We have studied this ques- 
tion carefully and have been able to show that there is 
definitely a T^ term in Kq/K arising from the infinite 
series of the ring diagrams, but we still do not know 
whether the T^ In T term of the exchange self-energy is 
exactly canceled by the dynamical effects arising from 
the infinite ring diagram series. Although the effect of 
the (T/Tp)"^ ln(T/Tp) term matters only at low temper- 
atures and low densities, where disorder effects dominate 
making it difficult to observe the log term experimen- 
tally, it is nevertheless important to establish whether 
this non-monotonicity associated with T^ In T survives 
higher-order correlation terms in the theory. We leave 
this as an unanswered theoretical question for the fu- 
ture. A direct experimental observation of the T^lnT 
term in the low-temperature inverse compressibility of 
2D electrons or holes will go a long way in settling this 
important question. 

To summarize, our goal of this paper is to understand 
the temperature dependence of compressibility in a high 
density homogeneous system where the interaction ef- 
fect is not too strong and HFA is valid. At very low 
density, disorder effects are of particular importance and 
the system may be highly inhomogeneous due to the for- 
mation of electron-hole puddles"*'^'^'®— liii^ii^Zr^, which 
arc not included in our theory. Our main result is that 
temperature effects in the compressibility could be quite 
important at higher temperatures. 
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Appendix A: Hartree-Fock compressibility in 3DEG 

In this Appendix, we provide the corresponding finite 
temperature Hartree-Fock compressibility for a standard 
three dimensional electron gas (3DEG). For 3DEG, the 

47re^ 

bare Coulomb interaction V{q) = 5-, where k is the 

background dielectric constant. We assume the total de- 
generacy g = 2 in the calculation. 

The finite temperature exchange energy of 3DEG is 
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jiven by 



Jn In - — 



e^x^-tLo) + l(Al) 



l/a;3 



e/3(efc/a;2-Ato) -|- 1 



where /3 = 1/fcsT, efc = fc^/2m, m is the effective 
mass, and /io is the finite temperature chemical potential 
without exchange energy. The non-interacting chemical 
potential, /io, is obtained by solving the following self- 
consistent equation, 



= ^ dx 



exp(- 

where /Jq = Ho/Ep, t — T/Tp and n — 

There is no explicit analytical formula for the finite tem- 
perature chemical potential in 3DEG. But at very low 
temperatures {t <^ 1), the asymptotic form of the non- 
interacting chemical potential is given by 



for 3DEG. 

37r2 



Mo « 1 - 
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(A3) 



The HFA chemical potential, /i, can be calculated by 
including the exchange energy for k = kp, i.e., 



fi = ^J.a(T) + T,ex{kp,T). 



(A4) 



Then, the normalized chemical potential fi = fi/Ep can 
be expressed 
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where a = (4/97r)^/"^rs with the dimensionless parameter 
Ts = (3/47rn)^/'^(me^/K), and 
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(A2) By differentiating Eq. (jA4l) with respect to n and using 



kp — iTT^n, the corresponding non-interacting compress- 
ibility of 3DEG is given by 
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(A8) 



where Kq and Kni are the zero temperature and the fi- 
nite temperature non-interacting compressibility, respec- 

tively, and /13 (t) — 1 -\ — . 

The finite temperature inverse compressibility Kq/K 
within HFA is given by 
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At zero temperature T 
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The HE chemical potential and the compressibility in 
the low temperature limit t <C 1 can be calculated by 
expanding the logarithmic function near x — 1. Then 
the chemical potential in Eq. (|A5|) becomes 

-(Aa+lnf) 
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(All) 



where A3 = 1 — 7 + In^'g) (where 7 is Euler's constant, 
with numerical value ~ 0.577216.). Finally we find that 



the low temperature 3D inverse compressibility within 
HFA as 



K 



2m2 
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(A12) 



The asymptotic behavior of high temperature compress- 
ibility is obtained by approximating the Fermi-Dirac dis- 
tribution function by the classical Boltzmann function. 
In the high temperature regime t — > 00, the normalized 
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chemical potential within HFA becomes 

H„,. = "°(57l;^)-|\/£-- 

and the corresponding high temperature inverse com- 
pressibility is calculated as 

which agrees with the non-interacting high-temperature 
classical result. 



The finite temperature chemical potential within static 
RPA, /i(T), can be calculated by including the self-energy 
of Eq. (HH), i.e., 

Ai(r) =/^o(T) + I]scif(fcf,r). (B5) 

At zero temperature (T = 0), the non-interacting part of 
chemical potential is the Fermi energy /i = Ep. For <C 
1, the chemical potential approaches the non-interacting 
value, which has been discussed in Sec.|TTl For interacting 
2D systems the self-energy within screened HFA has the 
asymptotic form (with gs = 2 and = 1) for ^ 1, 



Appendix B: Screened 2D Hartree-Fock 
compressibility 

We provide the screened Hartree-Fock compressibility 
(sometimes also called static RPA) for the 2D system, 
where the bare Coulomb interaction 'V^' appearing in the 
exchange self-energy is screened by the static RPA di- 
electric function. 

The self-energy with the finite temperature static RPA 
dielectric function is given by 



Sself(fc^^) = -y ^ 



£q V{q) 



<q,T) 



(Bl) 



where V{q) 



2^ 
Kg 



np (^q_fcy) is the Fermi distribution 



function. The static RPA dielectric function e {q, T) is 
given by^ 

e(q,T) = l + V-(g)H(g,r) = l + ^i^ (B2) 



where H((7. T) is the finite temperature static polariza- 
tion function and qs — 27re^H(g, r)/^ is the momentum 
dependent screening wave vector^. At low temperatures 
{T <t:Tp) the 2D polarizability becomes*. 



H(g,r)«^fl-e-^WT^, 



(B3) 



and its asymptotic form at high temperatures (T 3> Tp) 
becomes^ 
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Differentiating Eq. IB6I with respect to n and using the 
relation kp = 27rn, we get the zero temperature inverse 
compressibility within the screened HFA 



K{T = 0) 



(B7) 



Since is a positive value we have the positive com- 
pressibility even for 3> 1, which clearly disagrees with 
both the HE results and experimentsi"^. 

With Eqs. (jBl[) - (IB3|) . we find the asymptotic form of 
the self-energy at low temperatures (T ^ Tp) and for 
Ts > 1 
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We find that the leading order temperature dependent 
term in the self energy is exponentially suppressed. Con- 
sequently, the temperature dependent inverse compress- 
ibility at low temperatures becomes 
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At high temperatures (T ^ Tp), the inverse compress- 
ibility calculated within the static RPA approaches the 
HFA results given in Eq. ([TT]) because the non-interacting 
kinetic energy dominates. 
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